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ABSTRACT 

In order to try to understand the internal evolution of galaxies and relate this 
to the global evolution of the galaxy population, we present a comparative study 
of the dependence of star formation rates on the average surface mass densities 
(E M ) of galaxies at 0.5 < z < 0.9 and 0.04 < z < 0.08, using the zCOSMOS 
and SDSS surveys respectively. We derive star formation rates, stellar masses, 
and structural parameters in a consistent way for both samples, and apply them 
to samples that are complete down to the same stellar mass at both redshifts. 
We first show that the characteristic step-function dep endence of median s pecifi c 



star formation rate (SSFR) on Sjv/ in SDSS, seen by iBrinchmann et al.l (120041 ). 



is due to the changeover from predominantly disk galaxies to predominantly 
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spheroidal galaxies at the surface mass density logYiMchar ~ 8.5 at which the 
SSFR is seen to drop. Turning to zCOSMOS, we find a similar shape for the 
median SSFR - relation, but with median SSFR values that are about 5 — 6 
times higher than for SDSS, across the whole range of Ejv/, and in galaxies with 
both high and low Sersic indices. This emphasizes that galaxies of all types are 
contributing, proportionally, to the global increase in star formation rate density 
in the Universe back to these redshifts. The T, Mchar "step" shifts to slightly 
higher values of T, M in zCOSMOS relative to SDSS, but this can be explained 
by a modest differential evolution in the size-mass relations of disk and spheroid 
galaxies. For low Sersic index galaxies , there is little change in the size-mass 



relation, as seen by lBarden et al.l (120051 ) . although we suggest that this does not 



necessarily imply inside-out growth of disks, at least not in this redshift range. 
On the other hand, there is a modest evolution in the stellar mass - size relation 
for high Sersic index galaxies, with galaxies smaller by ~ 25% at z ~ 0.7. Taken 
together these produce a modest increase in T>Mchar- Low Sersic index galaxies 
have a SSFR that is almost independent of E^, and the same is probably also 
true of high Sersic index galaxies once obvious disk systems are excluded. 

Subject headings: galaxies: evolution, galaxies: high redshift 



1. Introduction 



One of the key unanswered questions in the study of galaxy evolution is what physical 
processes inside galaxies drive the changes in the star formation rates in individual galaxies 
that, taken together, prod uce the large de c line in the global star-formation rate density to 



redshifts since z ~ 2 (e.g. iLilly et al.lll996l ; iHippelein et al.ll2003l ; iHopkins fc Beacoml 12006 



and references therein) . It is known that there is a strong correlation between star formation 
rate (SFR) and stellar mass in the local Universe as shown by nume rous SDSS studies (e.g. 



Brinchmann et al. 2004: Salim et al 



2007 



Schiminovich et al.l 120071 ). This is demonstrated 



for example, in Fig. 17 of iBrinchmann et al.l (120041 ). although this figure also clearly shows 
that at high stellar masses the distribution of SFRs broadens significantly and the correlation 
between stellar mass and SFR breaks down. Many studies at intermediate redshifts of the 
SFR or speci fic SFR (SSFR) have been made as a function of the integra ted stellar mass of 
galaxies (e.g. iBauer et al.ll2005l : iBundy et al.ll2006l : iNoeske et al.ll2007al Jbl) but these did not 
use information on the internal structural properties of the galaxies. 



In contrast, studies using the local SDSS sample (e.g. IBrinchmann et al.l 120041 ) have 



argued that the surface mass density may be more important than stellar mass in regulating 
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star formation. Using the SDSS sample iBrinchmann et al.l (120041 ) found that the low specific 
star formation rate (SSFR) peak is more prominent at high Sjv/ than at high M*, and 
therefore concluded that the surface density of st ars is more impor t ant th an stellar mass 
in regulating star formation. In a follow-up study, iKauffmann et al.l (120061 ) found that the 
total spread in SSFR reaches a maximum at a characteristic stellar surface mass density 
^Mchar, and interpreted this as a qualitative change in the distribution of star formation 
histories above and below T, Mchar . The behaviour of the SSFR with E M follows a smoothed 
"step-function" dropping substantially at a characteristic Xa/c/w 

This observation and the resulting discussion therefore motivates the current investiga- 
tion. Using the HST/ACS images of the COSMOS field, plus star-formation rate informa- 
tion from emission lines measured in large numbers of zCOSMOS spectra we can study the 
changes that have occured in the SSFR-Sjv/ relation between redshifts approaching z ~ 1 
and the present epoch, as sampled by the SDSS studies, provided we can select comparable 
samples at the different redshifts, e.g. above a certain integrated stellar mass. 

The goal of this paper is therefore to obtain clues about the links between the internal 
evolution of galaxies, in particular the build-up of stellar mass, and the global changes that 
are seen in the population of galaxies as a whole, by studying this SSFR- Em relationship at 
significantly earlier epochs, 0.5 < z < 0.9. 

This work represents an improvement over an earlier study o f the S SFR versus stel- 
lar mass surface density in the COSMOS field by IZamojski et al.l (120071 ) by using secure 
spectroscopic redshifts, by deriving SFRs from emission line fluxes rather than ultraviolet 
continuum luminosities, by carefully constructing compatible mass-complete sub-samples for 
both zCOSMOS and SDSS, by computing half-light radii in a consistent way, at the same 
rest-frame wavelength, and by exploring the contributions of galaxies with different Sersic 
indices to the overall SSFR - relation. 

By focussing on internal properties of galaxies, this paper complements a number of 
other studies of the evolution of galaxies over the < z < 1 redshift range that are being car- 



ried ou t using the first 10,000 spectra from the zC OSMOS redshift survey (e.g. ICaputi et al. 



2008al lbl:l Mignoli et al.ll2008l ; ISilverman et al.ll2008l ). These other papers are more focussed on 
the integrated properties of galaxies and on the variation of these with the external galaxian 
environment. 

This paper is organized as follows: Sect. 2 describes our sample selection and the deriva- 
tion of the relevant physical quantities for both zCOSMOS (0.5 < z < 0.9) and SDSS 
(0.04 < z < 0.08) galaxies. In Sect. 3, we present our main results regarding the role of stel- 
lar mass and stellar mass surface density in regulating star formation activity for zCOSMOS 
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and SDSS galaxies as a function of morphology. These include new insights into the form 
of the SSFR - S^/ relation at the present epoch, and a comparison of the changes seen in 
this relation back to 0.5 < z < 0.9. These results are discussed in Sect. 3, and conclusions 
are presented in the last section. A cosmology with H = 70 km s -1 Mpc -1 , Q = 0.25, 
= 0.75 is used throughout this paper. 



2. Datasets 



2.1. The zCOSMOS sample at 0.5 < z < 0.9 



The Cosmic Evolution Survey (COSMOS. IScoville et al.ll2007al ) is the largest HST sur- 
vey ever undertaken, imaging an equatorial 1.64 deg 2 field with sin gle-orbit I-band exposures 
(IScoville et al.ll2007bl ). The zCOSMOS project (jLilly et al.l 120071 ) is securing spectroscopic 
redshifts for large numbers of galaxies in the COSMOS field. We now have in hand spectra 
of about 10 500 Iab < 22.5 I-band select ed galaxies over 1.5 deg 2 of the COSMOS field, the 
so-called zCOSMOS-bright 10k sample jLilly et alJliooi ). The zCOSMOS-bright sample is 
a flux-limited sample of galaxies with Iab < 22.5 generated from the COSMOS HST ACS 
images. The sample does not have a significant observational surface brightness selection, 
but, as with all samples selected on total magnitude, there is a selection in surface brightness 
at a given size. The success rate of measuring redshifts varies with redshift and is very high 
(more than 90%) between 0.5 < z < 0.9. For more details about the zCOSMOS survey, we 
refer the reader to iLillv et all ((2003, l2008h . 



Spectroscopic observations in zCOSMOS-bright were acquired using VIMOS with the 
R ~ 600 MR grism over a spectral range over 5550-9650A. This enables us to study rest-frame 
spectral features around 4000A, such as the [Oil] A 3727 emission line, to z ~ 1.2. 



2.1.1. Sample selection 



Only 10k sample galaxies with reliabl e redshifts a r e used for this study. Specifically, and 
with reference to the scheme decribed by ILillv et al.l (120081 ). we use redshifts of confidence 
classes 4, 3, 2.5, 2.4, 9.5, 9.4, 9.3, and 1.5. These have an overall reliability, based on repeat 
observations and on overall consistency with independent photometric redshifts, of more 
than 99%. 



We exclude about 500 stars, the broad line AGNs, and the few galaxies that lie outside 
the ACS fields, leaving us with a sample of 8131 galaxies at z < 1.2 for which structural 
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parameters are available 



in detail bv lSargent et al 



r om th e HST/ACS images from GIM2D Sersic fits, as described 
Automated morphological classifications using the ZEST 



scheme (ZEST; IScarlata et al 



20071 ) are also available. 



To have a high degree of spectroscopic completeness (see Fig. 9 in lLilly et al.ll2007l ) and 
for additional reasons discussed below, we further restrict the zCOSMOS sample for this 
study to galaxies lying between 0.517 < z < 0.900. 

To avoid night sky lines problems with our chosen tracer for star formation, the [Oil] 
emission line, we exclude all galaxies for which the center of the [O II] A 3727 line is less than 
12A away from the center of a strong night sky line. We therefore exclude galaxies with 
[O II] A 3727 lines close to the strong night sky lines 5893A, [OI]6300A, 6364A, or 6832A, i.e. 
at redshifts 0.5771 < z < 0.5852, 0.6871 < z < 0.6936, 0.7043 < z < 0.7108, and 0.8299 < 
z < 0.8363, respectively. We exclude galaxies in the redshift range z < 0.5174 to avoid 
problems with the [OH line flux measurements using the automatic software Platefit_VIMOS 
(ILamareille et al.ll2008l ) at the blue edge of the spectrum. This way we remain with 3232 
galaxies. 



2.1.2. Line measurements and AGN rejection 



Emission line fluxes were measured using the a utomatic routine Pl a tefit_V IMOS (ILamareille et al 



20081 ). After removing a stellar-component using iBruzual fc Charlotl (120031 ) models, Plate- 
fit_VIMOS completes a simultaneous gaussian fit of all emission lines using a gaussian profile. 
If an emission line measurement has a significance of less than 1.15<r (this value is derived 
using the distribution of fake detections, as described in ILamareille et al.ll2008l ). then an 
upper limit is calculated for the respective line. 

For galaxies for which Platefit_VIMOS detects all three emission lines [O III] A 5007, H/?, 
and [O II] A 3727, we use the [O III] A 5007/H/5 versus [O II] A 3727/H/5 diagram to distin- 
guish star formation dominated galaxies from objects obviously containing an active nucleus 
(narrow-line AGNs) wi th the aid of the empirical threshold derived using the 2dFGR S by 
Lamareille et al.l (120041 ) . Additionally, all the X-ray detected AGN ( iBrusa et al.ll2007l ) are 
excluded, regardless of their spectral properties. 

This results in a sample of 3048 galaxies at 0.5 < z < 0.9, comprising 1527 galaxies at 
0.5 < z < 0.7, and 1521 galaxies at 0.7 < z < 0.9. These objects are used for the analysis 
below. 
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2.1.3. Structural parameters 



Physical sizes are computed from the half-light radii, ri/2, defined to be the semi-major 
axis of the ell i pse c ontaining; half of the total flux, derived from the GIM2D Sersic fits of 
Sargent et al.l (120071 ). These GIM2D fits also provide Sersic indices. It should be noted that 
GIM2D provides a surface brightness profile for each galaxy as it would be in the absence of 
the instrumental PSF. 



As described in lSargent et al.l (120071 ). extensive tests and simulations showed that these 
fits slightly underestimate the half-light radii of galaxies. Therefore, a correction was applied 
to the half-light radii, following their Fig. 21, which depends on the half-light radius size in 
pixels, the I-band magnitude (i.e. the signal to noise), and the Sersic index of the galaxy in 
question. This correction was < 10% for 78% of the 3048 zCOSMOS galaxies at 0.5 < z < 
0.9, and more than 20% for just 2% of the sample. 



Although the ZEST automated morphological classification (IScarlata et al.l 120071 ) is 
available for the entire zCOSMOS sample, it has not been applied to the SDSS. There- 
fore, for consistency, we decided to simply use Sersic indices (n) to co mpare the structures 
(morphologies) of zC OSMOS and SDSS galaxies. As is well-known (e.g. lBlanton et al.ll2003al ; 
Scarlata et al.l 120071 ) n > 2.5 galaxies correspond mostly to galaxies classified as early-types 
spheroid-dominated galaxies, while n < 1.5 galaxies are mostly objects classified as late-type 
disk-dominated galaxies. 

In what follows, we will split the samples according to Sersic index, and loosely re- 
fer to these as disk- and spheroid-dominated galaxies, as late and early type etc., while 
endeavouring to keep the Sersic criteria clear. 



2.1.4. Ellipticity 



In addition to the Sersic indices, we also use an ellipticity (q 
separate galaxies. There are two motivations for this: 



1 — b/a) criterion to 



First, iMaller et al.l (120081 ) have shown that 90% of galaxies with an axis ratio b/a < 
0.55 are disk galaxies in an SDSS sample. Therefore, eliminating high Sersic galaxies that 
nevertheless have a low b/a ratio should produce a cleaner sample of early-type galaxies, 
since these highly inclined galaxies are unlikely to be true early types. Secondly, the low 
Sersic ind ex disk galaxies with b /a < 0.55 are likely to be more affected by dust extinction 
(see, e.g.. iMollenhoff et al.ll2006l ). affecting estimates of stellar masses based on the spectral 
energy distributions, star formation rates derived from emission lines and possibly also their 
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size estimates, all of which become more difficult. Eliminating highly inclined galaxies with 
low Sersic indices also therefore has advantages. 

For both these reasons, we will often focus on the subset of galaxies with b/a > 0.55 in 
both the zCOSMOS and SDSS samples, and will do so regardless of their Sersic index. By 
doing so, we obtain a cleaner sample of early types and eliminate the hardest to interpret 
late types. Applied consistently between the SDSS and high redshift samples, neither should 
introduce systematic biasses. 



2.1.5. Stellar masses and surface mass density 



Stellar masses, and the surface mass densities, must be derived from estimating the 
mass-to-light ratio of the stellar population using multi-band photometry. The advantages 
of doing this using the widest possible wavelength range and extending as far as possible 
into the rest-frame near-infrared are well-known. 

Unfortunately, while such data exist in the COSMOS field, systematic near-infrared 
photometry is not yet available for SDSS. Therefore, in order to achieve the highest possible 
internal consistency for our study, we therefore derive mass estimates from purely rest-frame 
optical colors, for both SDSS and zCOSMOS, using the equation given below in Eq.[TJ 

We have however compared thes e stellar masses with mas ses derived using the entire 



COSM OS optical to infrared SED and lBruzual fc Charlotl (120031 ) models by lBolzonella et al. 



( 120081 ) . As shown in Fig.[T] these show good agreement with a statistical rms of around 



as described by 


Bolzonella et al. 


Maraston 


2005) would produce 



a systmatic shift of ~ 0.1 dex towards lower SED masses than those shown in Fig.[H thus 
eliminating the offset between the masses computed using Eq.UJand the SED masses. 

Stellar masses are derived using the relation b etween res t -fram e U-B and B-V colors and 
mass-to-light ratio (M/L), using the Eq. (1) from lLin et al.l (120071 ). which corrects the M/L 
for e yolution, and acc ounts for variations in the mass-to-light ratio with color. The relation 
from lLin et al.l (120071 ) is modified as follows, and applied to both SDSS and zCOSMOS data 
sets: 
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log(M*/M e ) = 0.4 • (5.48 — Mb — 0.10) + 1.737 • (M B - M v + 0.1) 
+0.098 • {M v - M B - 0.81) - 0.130 • (M v - M B - 0.81) 2 

-0.268 • z - 1.003 + log(1.7) 



(1) 



The difference between Eq. (1) in lLin et al.l (120071 ). which uses Vega magnitudes, and our 
Eq. Uis due to th e trans formation onto the AB system, using the con versions given in Table 1 
in lWillmer et al.l ( 20061 ) . Fin ally, since our e stimate of the SFR from lMoustakas et al.l ( 2006) 



is based on a Salpeter IMF JSalpjtejjl]jp) . and the stellar masses computed by iLin et al.l 



( 120071 ) assume a Chabrier IMF (jChabrierl 120031 ) . we convert the latter m asses to a Salpeter 
IMF statistically, using the conversion factor of 1.7 (IPozzetti et al.l 120071 ). This accounts for 
the last term in Eq.[TJ 

For zCOSMOS galaxies, M Uy M B , and My are Johnson rest-frame absolute magnitudes 
in the AB system and were calculated from the best-fitting continuum template for each 
galaxy. This was obtained by fitting the multi- wavelength COSMOS spectral energy dis - 



tribution (SED) that is available for each object (Capak et al.l 120071 ; iTaniguchi et al.l 120071 ) . 



using the ZEBRA photo-z code (IFeldmann et al. 



200 



(]) at the known, spectroscopic, red 



shift. To achieve internal consistency with the size measurements, the SEDs were previously 
normalized to the actual (half-light) /-band magnitude measured from the GIM2D fits. 

Finally, surface mass densities were computed using logS M / '[M Q kpc~ 2 ] = log(M*/ M Q ) — 
log[2ir(ri/2/ 'kpc) 2 ], for both the zCOSMOS and SDSS galaxies. Further details on the cal- 
culation of half light radii for the SDSS sample are given in Sect. 12. 21 Since many studies 
of galaxy evolution at z < 1 are based on stellar masses, in Fig.[2]we show, for the benefit 
of the reader, the comparison between Em and stellar mass for the zCOSMOS and SDSS 
mass-complete samples (see below the discussion of the mass completness) . A slight increase 
of T>m with stellar mass, but with a large scatter, can be seen for both zCOSMOS and SDSS 
samples of relatively massive galaxies. 



2.1.6. Star formation rates 



One of the most reliable and well-calibrated star formation rate (SFR) indicators is the 
Ha emission line, but the drawback for many high redshift studies is that Ha is redshifted 
out of the optical window beyond z ~ 0.5, and its measurement req uires very time-consumin g 



near-infrared spectroscopy for small samples of galaxies (see e.g. iMaier et al.l 120051 . 120061 ). 



Therefore, the strongest emission feature in the blue, the [O II] A 3727 forbidden-line doublet, 
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is usually used at z > 0.5 as a tracer of the SFR. We use the [O II] A 3727 to SFR calibration 
of lMoustakas et all (bood ) for this study, for both the zCOSMOS objects at 0.5 < z < 0.9 
and the SDSS galaxies at 0.04 < z < 0.08. The two dominant sources of scatter in optical 
SFR diag nostics, dust extinction and metallicity, both correlate strongly with L(B) (see e.g. 
Fig. 16 in lMoustakas et al.ll2006l ). and this motivated their empirical calibration of SFR in 
terms of the [O II] A 3727 luminosity. This calibration was developed by tuning [OII]-derived 
star formation rates to match those based on extinction-corrected Ha using observations of 
nearby galaxies. It might be a concern that the extinction changes at intermediate redshift, 
but our own obser vations of Ha in a limited sample of CFRS galaxies at 0.5 < z < 0.9 
( IMaier et al.l 120051 ) indicate that neither average extinctio n nor average metallic ity have 
greatly changed for most galaxies at this redshift. Moreover, iMoustakas et al.l (120061 ) showed 
their calibration to be effective also at intermediate redshifts. 

In detail, the conversion of [O II] A 3727 luminosity to SFR is undertaken as follows. 
First, because the slits in the VIMOS masks have a width of 1 arcsec, an aperture correction 
to the [O II] A 3727 line flux is applied before computing the SFRs. Each zCOSMOS spectrum 
is convolved with the ACS 1(814) fil ter and then this ma gnitude is compared with the I- 
band magnitude of the GIM2D fits of ISargent et al.l (120071 ). The difference between the two 
magnitudes gives the aperture correction factor for each spectrum. For more than 85% of 
the sample, these aperture correction factors are between 1 and 3. This correction assumes 
that the [O II] A 3727 line flux and I-band continuum suffer equal slit losses. 

Then, the corrected [Oil] line luminosities are transformed into SFRs using a correction 
factor based on the galaxy's B-ba nd absolute magn i tude, as given by a linear interpolation 
of the < R > values in Table 2 of lMoustakas et al.l (120061 ). and shown in their Fig. 19. We 
use the followin g formula to convert the [Oil] luminosity into SFR assuming a Salpeter IMF 
( ]Salpeterlll955l ) for both luminous SDSS and zCOSMOS galaxies: 



logtSFR/Meyr" 1 ] = log [L [0 n]/ (ergs s" 1 )] - 41 - 0.195 ■ M B - 3.434 



(2) 



We h ave checked tha t the S FRs derived using Eq.[5] are in broad agreement with those 
derived by lPozzetti et al.l (120081 ) using the entire COSMOS optical to infrared (including K, 
3.6/xm, and 4.6/im) SED. 



2.1.7. Mass completeness of the zCOSMOS sample 



Both zCOSMOS and SDSS are flux-limited samples. Since the mass-to-light ratio will 
depend on the stellar population, which itself will produce different SEDs and rest-frame col- 
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ors, both the SDSS and zCOSMOS samples will have a complicated mass-selection function. 
The problem is more acute for zCOSMOS because we can choose the SDSS redshift range so 
that the SDSS sample is effectively complete for the masses of interest. Of course, working in 
terms of stellar mass instead of luminosity is useful not only because it is closer to quantities 
predicted by theory, but also because it in principle removes the evolution in luminosity due 
to the aging of stellar populations. Blue late-type star-forming galaxies generally have lower 
mass-to-light ratios. Therefore, in a given flux-limited survey such as zCOSMOS, they are 
detected to smaller masses than are red early-type galaxies, simply because the latter are 
fainter at a given mass, i.e. have a higher mass-to-light ratio. This is shown in Fig.[3] (stellar 
mass versus redshift) and FigJH (SSFR versus stellar mass diagram). 

For this reason, if we calculate the maximum mass-to-light ratio at a given redshift, i.e. 
that of a red passively evolving galaxy, and apply this to calculate the minimum mass that 
a galaxy must have in order to be brighter than our flux limit, then we can safely assume 
that all other galaxies with masses above this mass limit will be visible within the survey, 
because they will have lower mass-to-light ratios. 



We have explored this by running a wide range of iBruzual fc Charlotl (120031 ) models, 
with different SFR e-folding timescales (r from 5 Gyrs to 0.7 Gyrs), and with different starting 
redshifts (zj from 1.3 to 5). For each model, we scale the model so that it has an observed 
Iab = 22.5 (the flux density limit of zCOSMOS) at either z = 0.7 or z = 0.9, and plot the 
mass and SFR of this "just-detectable" model on Fig. @] as large magenta dots. While there 
is scatter because of the range of models used, these magenta dots give a good indication of 
the mass completeness limit of the sample: galaxies that occupy the diagram to the left of 
the filled circles will not be seen while those to the right will have been detected. 

As the mass-to-light ratio asymptotes to a limiting value the locus of the dots becomes 
vertical and this limit is adopted as the overall mass completeness limit. These adopted 
mass completeness limits are shown by the green solid vertical thick lines in the two panels 
of Fig.H 

Henceforth, we only consider objects above these mass limits for the two redshift zCOS- 
MOS redshift ranges. At logM^ > 10.4, we obtain a mass-complete zCOSMOS sample of 
648 galaxies at 0.5 < z < 0.7, while for logM^ > 10.7 we obtain a mass-complete zCOSMOS 
sample of 520 galaxies at 0.7 < z < 0.9. As noted above and described below, the SDSS 
redshift range is chosen so that the SDSS objects are intrinsically faint enough that we can 
obtain the equivalent mass-complete samples also in SDSS. In what follows, we will always 
compare zCOSMOS and SDSS samples selected to the same mass limits. It should be noted 
that these mass limits are determined at the epoch at which the galaxy is observed, and do 
not, therefore account for any stellar mass added in between. As we remark below, we would 
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expect this added mass to be modest (0.1 dex) in most cases. 



2.2. The SDSS comparison sample 

We have selected a local comparison sample of SDSS galaxies from the DR4 release, 



Garching repository ( http: / / www.mpa-garching.mpg.de/SDSS/DR4 / ) , in the redshift range 



0. 04 < z < 0.08. The l ower redshift limit of z — 0.04 is chosen following the recommendation 



of iKewley et al.l (120051 ) to reduce the systematic and random errors in the SFRs that arise 
from aperture effects due to the 3arcsec size of the SDSS fibers. The upper limit choice of 
z = 0.08 ensures mass completeness of red (high mass-to-light) SDSS galaxies at logM^ > 
10.4, as shown in panel b) in Fig.[3j We selected galaxies with Petrosian r magnitudes in the 
range 14.5 < r < 17.77. The bright limit is necessary because SDSS becomes incomplete for 
bright galaxies of large angular size and to avoid objects with saturated SDSS photometry, 
whereas the faint limit corresponds to the nominal magnitude limit of the main galaxy sample 
in SDSS. It should be noted that duplicate objects, and SDSS galaxies on problematic or 
special plates are excluded from the SDSS sample. This exclusion is based solely on position 
on the sky and should not be related to any galaxy properties. 

To ensure a consistent comparison of the physical properties of the SDSS and zCOSMOS 
samples, we perform the following steps to derive stellar masses, sizes, Sersic indices, and 
SFRs for SDSS galaxies: 

a) To derive the sizes r e i/ 2g of SDSS galaxies, we start from the Petrosian radius con- 
taining 50% of the Petrosian flux in the g-band, rpet 50g . The observed g-band for SDSS 
galaxies at 0.04 < z < 0.08 corresponds approximately to rest-frame B-band (to be precise, 
effective wavelengths between 4587-4417 A) and is therefore consistent with the derivation of 
sizes for the zCOSMOS sample from the observed ACS I-band images, which have effective 
wavelengths of 5373A at z = 0.5, 474lA at z = 0.7, and 4242A at z = 0.9. Since these are 
rather similar, and always above the 4000A break, any dependence on redshift (the so-called 
"morphological k -correction" ) s h ould be small. We transform rpet 50g into half-light radii 



using Eq. (6) of I Graham et al.l (120051 ): ri/2 9 = rpet^o g /(l — 6 x 10 6 (rpet$ 0g /rpet 5 - 



b 50g) 

where rpetgo g is the Petrosian radius containing 90% of the Petrosian flux. This equation 
uses the concentration rpet 90g /rpet 50g to correct for the light outside the Petrosian aperture: 
this flux deficit is 0.20 mag in the case of de Vaucouleurs R 1 ^ profile, and 0.50 mag for an 
i? 1 / 8 profile. Then, to be consistent with the zCOSMOS r\/ 2 measurements and since the 
SDSS Petrosian flux was measured with a circular aperture, we transform r\j 2g to r e i/ 2g , 
the semi-major axis of the ellipse containing half of the total flux: r e i/ 2g = T\i 2g x yja/b, 
where a/b is the ratio between the major and minor axis of the ellipse. We did several tests 
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to assure the consistency of the derivation of sizes for the zCOSMOS and SDSS samples, 
and the results are presented in appendix |A] This r el / 2g value, the SDSS half-light radius 
shown in the following in different diagrams, is then used to derive surface mass densities as 
described in Sect . 12. 1.51 



b) For the reason explained above in Sect. 12. 1.31 we use Eq.[T] to derive stellar masses 
for the SDSS galaxi es. SDSS u, g , andr (Petrosian) magnitudes have been k-corrected to 
z = according to iBlanton et al.l (j2003bl ). using k-correct version 4.1.4 which is available 
at http:/ /cosmo. nyu.edu/blanton/kcorrect/ U-, B- and V-band rest-frame AB magnitudes 



w ere computed from these k -corrected u, g and r magnitudes, using the relations from Table 2 



m 



Blanton &: Roweid (120071 ). and stellar masses were calculated using these values in Eq.[TJ 



The stellar masses obtained in this way were found to be i n good agreement with the SDSS 
total stellar masses estimated from the SDSS spectra by iKauffmann et al.l (j2003al ). The 
comparison shows a statistical rms of around 0.16 dex per galaxy and an offset of 0.03 dex 
i n the mean. This is com parable to the quoted uncertainties in the spectral measurements 
(IKauffmann et al.ll2003al), and with the scatt er and systematic differences between different 
spectral mass estimates ( Gallazzi et al. 20051 ). 



c) We calculate SDSS SFRs from the [O II] A 3727 emission line fluxes using Eq.EJ 
the same equation used for zCOSMOS galaxie s. These SDSS SFRs are then additionally 
corrected using the aperture corrections given bv lBrinchmann et al.l ( 20041) . corrections which 
should work for z > 0.04, as also demonstrated by iKewley et al. "(|2005h . 



d) We use the Sersic indices in the g-band from the New York Univ e rsity value-added 
galaxy catalog (NYU-VAGC), from Sersic fits described in IBlanton et all tood ). 



e) We exclude AGNs using the [O HI] A 500 7/ Hff versus [N Hi A 6584 /Ha diagnostic 



diagram, excluding objects that satisfy Eq. (1) of Kauffmann et al. (I2003bl ). 



This results in a sample of 55230 SDSS galaxies at 0.04 < z < 0.08. Applying the same 
mass cut as applied to produce the two mass-complete zCOSMOS samples, we are left with 
21497 SDSS galaxies at 0.04 < z < 0.08 with logM, > 10.4, and 8904 SDSS galaxies at 
0.04 < z < 0.08 with logM, > 10.7. 



3. Results 



3.1. Specific star formation rates and downsizing 



The star formation rate per unit stellar mass, the specific star formation rate, is an 
indicator of the galaxy star formation history, since 1/SSFR defines a characteristic timescale 
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of the stellar mass build-up. The inverse of the specific SFR, Tsfr, is the time required for 
the galaxy to form all of its stellar mass at the current SFR and this is shown on the right 
axis of Fig.|H The ages of the universe for the two respective redshift ranges considered 
are shown as the horizontal dashed magenta lines. Galaxies with Tsfr higher than the age 
of the Universe at that time are relatively "quiescent" and must have had a higher star 
formation rate in the past, while those with Tsfr less than the age of the Universe may 
be thought of as "active" or "forming" , in the sense that they cannot have maintained this 
high star formation rate for all of their lifetime. If the star formation rate is constant or 
slowly declining, then "forming" galaxies will follow a roughly diagonal track in FigJUas they 
age towards lower redshifts, whereas "quiescent" galaxies will follow a nearly vertical path 
downwards towards lower SSFRs (the two models on the right in the panels of FigHJ). 

It should be remembered that Fig.H] shows all galaxies in the zCOSMOS samples, re- 
gardless of whether they belong to the mass-complete samples to the right of the heavy 
vertical lines in the two panels. This figure shows two well-known phenomena: first, the 
well-known downsizing effect is observed in the sense that at 0.5 < z < 0.7 very few galax- 
ies with masses above 1O 1O ' 8 M (to the right of the vertical, magenta line, in the yellow, 
hatched region) have SSFRs above the magenta lines, while at 0.7 < z < 0.9 there exist 
several dozens of galaxies above 1O 1O,8 M that are "forming". Secondly, it is also obvious 
from the figure that galaxies with different Sersic indices have different star formation his- 
tories, even at a fixed mass: late-type galaxies with n < 2.5 have mainly high SSFRs, which 
indicate continuing star formation, and possibly a recent onset of star formation, whereas 
early-type galaxies with n > 2.5 have mostly low SSFR with Tsfr higher than the age of 
the Universe, indicating higher SFRs in the past and/or a higher redshift of their onset of 
star formation. 



3.2. The SDSS SSFR - E M relation 

In Fig.0, we plot the SSFR versus T, M diagram for the 0.04 < z < 0.08 SDSS galaxies 
with logM„ > 10.4 - the sample that should be directly comparable to the 0.5 < z < 0.7 
zCOSMOS sample. Individual measurements are shown as green dots, t he median SSFR 



values in different Em bins as magenta triangles, and, following Fig. 9 of iKauffmann et al. 



( 120061 ). 25th and 95th percentiles of the distribution of SSFR in the respective Y> M bin as 
solid magenta lines. The error bars shown are estimates of the error in the median, and are 
computed as Ax = (xo.84 — xo.xe) / 'VjV : , where N is the number of galaxies in each E M bin, 
and x .84 and £ .i6 denote the 84th and 16th percentiles of the SSFR distribution. Panel a) 
shows all galaxies, while panel b) shows only the SDSS galaxies with an axis ratio b/a > 0.55 
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(see 12.1.41 above). Panel c) shows the low Sersic index (n < 1.5) objects with b/a > 0.55, 
and panel d) shows the n > 2.5 SDSS galaxies, again with b/a > 0.55. 



In a previous study of the SSFR - Em relation for the SDSS sample, iKauffmann et al. 



( 120061 ) found evidence that the distribution of star formation histories changes qualitatively 



above a chara cteristic surface mass d ensity logY^Mchar ~ 8.5. Moreover, for higher mass 



SDSS galaxies. IKauffmann et al.l (120061 ) showed in their Fig. 9 that the average SSFR remains 
constant below T> Mchar and decreases at higher surface mass densities. This behaviour is seen 
in panels a) and b) of Fig.[5j 

This trend is clearly due to the change-over of different structural types from disk- 
dominated low Sersic index (n < 1.5) galaxies (in panel c of Fig. [5]) to bulge-dominated high 
Sersic index (n > 2.5) galaxies (panel d of Fig.E]) as the E M increases. This change-over is 
clearly seen in Fig.JHl which shows the fraction of n < 1.5 and n > 2.5 SDSS galaxies (always 
with b/a > 0.55) in different Em bins (magenta open and filled triangles). The fraction of 
n > 2.5 SDSS objects shows a sharp increase at the point where the SSFR abruptly changes 
in Fig.0 i.e., at a transition surface mass density /o(7S A / irans ~ 8.45 (dashed magenta vertical 
line in Fig.E]), close to the characteristic surface mass density logY, Mchar ~ 8.5 of the break 
in the SSFR - Em relation. 

Inspection of panel c) of Fig. [5] shows that the median SSFR for the late-type low Sersic 
SDSS galaxies {n < 1.5) and the plotted percentiles are all almost constant with Em - they 
decline very slightly with increasing Em- In contrast, we see a rather steep drop in both the 
median and the 25th percentile of SSFR for the early-type n > 2.5 SDSS galaxies (panel d), 
down to a SSFR level about 5-6 times below that of the late-type galaxies, where the bulk 
of the n > 2.5 sample resides. We return to the form of the SSFR— Em for n > 2.5 galaxies 
(panel d) below. In the meantime, we note that Fig.Oand panels c) and d) of Fig. explain 
how the shape of the median SDSS SSFR— Em relation for all galaxies is the result of the 
different SSFR— Em relations for n < 1.5 and n > 2.5 galaxies and the increasing fraction 
of early- type n > 2.5 SDSS galaxies with increasing Em above T^Mtrans- The large difference 
in the median SSFR for early and late-type (n > 2.5 and n < 1.5) galaxies "automatically" 
makes T> Mchar ~ T> Mtrans . 

This explanation of the overall SSFR— Em relation should no t be surprising but has not , 



as far as we are aware, been remarked upon previously, although ISchiminovich et al.l (120071 ) 



did note the different behavior of n < 2.5 and n > 2.5 SDSS galaxies. 
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3.3. The zCOSMOS SSFR - E M relation at 0.5 < z < 0.9 

Figures [7] and [8] show the SSFR vs. E^ for zCOSMOS galaxies at 0.5 < z < 0.9, and 
the comparison with the median and percentile curves summarizing the SDSS relations from 
Fig. [51 The individual measurements for the zCOSMOS mass-complete sample at 0.5 < z < 
0.7 (Fig.EJ) and 0.7 < z < 0.9 (Fig.ED are shown as black dots, and zCOSMOS median SSFR 
values for different Em bins as cyan, filled squares. 

It should be noted that in SDSS the number of galaxies with no measurements (upper 
limits) of SSFR is always lower than 50% in all bins we are considering, so the median 
SSFR in the respective Em bin is always well-defined. Unfortunately, this is not the case for 
the highest Em bins (logT^M > 9) in the zCOSMOS sample: more than 50% of the galaxies 
in these bins have only upper limits for the [Oil] fluxes from Platefit_VIMOS. To overcome 
this problem, we first simply computed a median by setting all upper limits to actual values. 
However, we also generated an "average" co-added spectrum of all the spectra in each bin, 
and derived a mean SSFR from that, as follows. Each individual zCOSMOS spectrum (with 
a spectral resolution of 2.55A) was first transformed to its rest-frame and rebinned to a 
common grid of 2 A per pixel. These individual rest-frame spectra were then normalized 
to the average value in a featureless region of the continuum and averaged. The [Oil] flux 
from the co-added spectrum was transformed to a SFR using the usual calibration, and to 
a SSFR using an average mass of the galaxies in the respective E^ bin. In fact, the mean 
SSFR value for the logT*M > 9 bins (shown as green open squares in Figs. [7] and [8]) derived 
in this way were found to be in good agreement with the median values that were obtained 
by setting the upper limits to their maximum values (cyan symbols). 

For reasons mentioned in Sect . 12.1."^} we focus primarily on galaxies with b/a > 0.55. 
Comparison of panels a and b of Figs. [7] and [8] shows that this does not have a big effect 
on the analysis. The SSFR— Em relation in zCOSMOS has the same general shape as that 
for the SDSS across the whole range of Em- However, the curve is shifted across the full 
range of Em to SSFR values that are about a factor of 5—6 higher in zCOSMOS than in 
SDSS. There is also a small shift of about 0.1-0.2 dex in the characteristic T^Mchar towards 
higher surface mass densities at higher redshifts, indicated by a localised blip in the change 
in SSFR around T^Mchar- The shift in T^Mchar is seen in the small panels of Fig. [7] which show 
the difference between the median zCOSMOS and SDSS SSFR values in a given E^ bin. If 
the shape of the SSFR-Em relation stays the same but there is a horizontal shift in Emc/md 
then the difference between the median zCOSMOS and SDSS SSFRs will be constant but 
with a localized excess in the region over which E Mchar has increased. This is seen as a blip 
in the small panels a, b, and d of Fig.[7J To illustrate more clearly the amount by which 
^Mchar shifts (0.1-0.2 dex) towards higher surface mass densities at higher redshifts, we have 
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also included in Fig. 7 A an insert showing a part of panel a with a finer grid in Em- 

The rise in median SSFR is seen both in the entire sample of zCOSMOS and SDSS 
galaxies (with or without the b/a > 0.55 cut) and in the individual n < 1.5 and n > 2.5 
subsets shown in panels c and d of the Figs. [7] and [HJ This rise in SSFR is essentially the 
same as the rise in th e global star-formation rate density (SFRD) in the Universe as a w hole 



to this redshift (e.g. iLilly et al.lll996l ; iHippelein et al.ll2003l ; iHopkins fc Beacoml 120061 . and 
references therein). The importance of this result is that it indicates that galaxies of all 
are contributing, proportionally, to this global evolution in the SFRD. Although the high 
Em galaxies with Em > SmcW (mostly early type galaxies with n > 2.5) have a SSFR that 
is an factor of about six lower, on average, than the galaxies of the same mass with lower 
Em < ^Mchar (generally late type galaxies with n < 1.5), the rise in SSFR by a factor of 5-6 
is seen across the board. 

As in SDSS, the SSFR for n < 1.5 galaxies (panels c in Figs. [7] and [8]) remains roughly 
constant with Em, albeit at a level that is elevated with respect to SDSS by about a factor 
of six. The shape of the relation for n > 2.5 galaxies (panel d of Figs. [7] and [S]) is also similar 
to that seen in the SDSS, with a decline in the SSFR around T^Mchar and a roughly constant 
SSFR at higher surface mass densities. The SSFR is again uniformly elevated across the 
range of Em by a factor of about five to six. 

With the better resolution of the ACS COSMOS images we can now examine the galaxies 
which are causing the upturn in SSFR of n > 2.5 galaxies at the lower Em in panel d of 
the Figs. [7] and [HJ Most of the zCOSMOS galaxies with high SSFRs and lower Em that are 
responsible for this upturn are clearly disk galaxies with a dominant bulge. If we remove 
these from the sample, as indicate by the red points in the panel d of Fig.[7J then we find 
a SSFR— Em relation that is essentially flat. We suspect that the same would also hold for 
the n > 2.5 SDSS sample. Although we have had to parameterize the galaxy morphology in 
terms of a crude Sersic index, in order to achieve uniformity between SDSS and zCOSMOS, 
we would predict that the form of the SSFR— Em would consist of two essentially flat relations 
in SSFR— Em for disks and spheroids separately, which partly overlap in Em and which are 
offset in SSFR by a factor of five to six. Both components appear to increase in SSFR by 
the same amount to the redshifts probed by this study. 



3.4. The SDSS and zCOSMOS stellar mass-size relations 



In order to better understand the small shift in Y>Mchar with redshift, we have looked at 
the relationship between mass and half-light radius in both Sersic subsamples at both SDSS 
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and zCOSMOS redshifts. Figure [9] shows the semi-major half-light radius in the observed 
ACS I-band (corresponding roughly to rest-frame B-band for this redshift range) versus 
stellar mass for the zCOSMOS mass-complete sample at 0.5 < z < 0.7. Black dots are 
individual zCOSMOS measurements, while cyan squares denote the median ri/2 values in 
the respective mass bin. The SDSS median half-light radii in different mass bins derived 
from g-band images (as described in Sect. 12. 2p are shown as filled, magenta triangles. The 
16th and 84th percentiles of r\i2 i n the SDSS and zCOSMOS galaxies in each mass bin 
are shown as solid magenta and cyan lines respectively. The two dashed black lines show 
constant mass surface densities of logY^u = 9 (upper line), and logY^u = 9.5, respectively. As 
before, panel a) of Fig. [9] shows all galaxies in the zCOSMOS mass sample, panel b) galaxies 
with b/a > 0.55, panel c) galaxies with b/a > 0.55 and n < 1.5, and panel d) galaxies with 
b/a > 0.55 and n > 2.5. The distribution of sizes of SDSS galaxies is displaced to larger sizes 
at a given mass compared with zCOSMOS when galaxies of all Sersic indices are considered 
(panels a and b). However, inspection of panels c and d of Fig. [9] shows that this effect is 
almost entirely due to galaxies with n > 2.5. The respective stellar mass-size relation for 
the slightly higher redshift 0.7 < z < 0.9 galaxies is shown in Fig.[10j and is discussed at the 
end of this section. 

Panel c of Fig.[9] shows that the stellar mass-size relation of n < 1.5 zCOSMOS galaxies 
at 0.5 < z < 0.7 does not differ significantly from t he similar rela t ion fo r SDSS n < 1.5 



galaxies. This is consistent with the GEMS result of iBarden et al.l (120051 ). These authors 
interpreted the almost negligible evolution with time in the stellar mass - size relation as 
support for an "inside-out" scenario. We would interpret an "inside-out" scenario to involve a 
broadly constant surface mass density (whether measured at the center or at a characteristic 
half mass radius) coupled with an increasing half-mass radius, i.e. an evolution parallel to 
the dashed lines in Fig. [91 

However, as remarked above, the SFR in these galaxies declines rapidly (by a factor of 
about six as shown in panel c of Fig. [7]) since z ~ 0.7. The subsequent evolution of these 
n < 1.5 zCOSMOS galaxies is described presumably by a declining SFR, as suggested by 
the paths of model galaxies shown in Fig.H] (solid green tracks). We show in panel c of 
Fig. [7] the possible tracks of model galaxies (red solid arrows) with declining SFRs and with 
sizes assumed to be constant over this time period. The mass, and thus stellar mass density 
increases by AM* < 0.3M*. Given that the stellar mass-size relation seen in panel c of Fig. [9] 
is almost flat, it is difficult to see much evidence for the required diagonal evolution in the 
size-mass plane. We therefore do not believe, over this redshift interval at least, that there 
is much evidence for inside out growt h of disks. Thi s is also consistent with the lack of 



evolution in the size function of disks (ILilly et al.lll998t ISargent et al.ll2007l ) 
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In contrast, panel d in Fig. [9] illustrates the difference in the stellar mass- size relation 
of high Sersic index (n > 2.5) zCOSMOS galaxies compared with the corresponding SDSS 
sample. Since the masses of early-type n > 2.5 galaxies are not expected to evolve signif- 
icantly between z ~ 0.7 and z ~ because of the low and rapidly declining SSFRs (the 
two models on the right in the left panel of FigfjJ, we conclude that zCOSMOS objects are 
on average 25% smaller at a given mass (panel d and D in Fig. [9]). This change in size is 
similar to but smaller than found by lTrujillo et all (120071 ). who claimed ~ 60% smaller ri/ 2 
on average compared to SDSS for objects at similar redshift with somewhat higher stellar 
masses (logM* > 11). 

One might worry that stellar masses for n > 2.5 zCOSMOS galaxies using Eq.[T] are 
systematically wrong. To check the reliability of the calculated zCOSMOS masses, we com- 
pare the B-band mass-to-light ratios of the zCOSMOS n > 2.5 galaxies with the M/L ratios 



obtain ed by studying the fundamental plane for early-type (E+SO) galaxies by iTreu et al. 



( 120051 ). Figure [TT] shows the B-band M/L ratios as a function of mass for our zCOSMOS 
sample (black f illed circles ), the M/L ratios of spheroids (E+S0) galaxies at 0.5 < z < 0.7 
from Fig. 15 of ITreu et al.l (120051 ) . and the local relation for early-type (spheroid) galaxies 
(dashed line) taken from their same figure. The slope of the mean observed zCOSMOS M/L- 
mass relation agrees with the slope of the local relation, as indicated by the solid diagonal 
line in Fig. [Til which is obtained by shifting the local relation to the lower M/ L expected 



for a p assively evolving population. Moreover, the 0.5 < z < 0.7 spheroids from ITreu et al. 



( 120051 ) occupy a similar region of the M/L - mass diagram as the zCOSMOS galaxies at the 
same redshifts, suggesting that our stellar masses derived from optical colors for n > 2.5 
(early-type) galaxies are reasonable, and further suggesting that the changes in the stellar 
mass-size relation at 0.5 < z < 0.7 in Fig. [9] are due mainly to smaller sizes (or at least 
smaller measured 7*1/2), a t a given mass, in zCOSMOS compared with SDSS. 

Similar r e sults regarding the size evolution of early-type galaxies have been seen by 



Truiillo et al.l (120071) . with an even larger change in size, and at higher redshifts z > 1.4 by 



Cimatti et al.l (120081 ). Both these authors advocated dry mergers as a p ossible mechanism 



for the growth in size. Simulations with realistic boundary conditions by lNaab et al.l (120071 ) 
indicated that mergers could increase the size of galaxies whi le the m ass changes only slightly. 
In the smooth envelope accretion scenario of iNaab et all (120071 ). accreted stars (mainly 
provided by min o r mer gers) form an envelope whose size increseases smoothly with time. 



Khochfar fc Silk! (120061 ) used a semi-analytical model of galaxy formation to predict the 



redshift-size evolution of elliptical galaxies, and their Fig. 4 shows that elliptical galaxies 
with masses above 10 10 M o could increase their sizes by a factor of about 1.25 from z ~ 0.7 
to z ~ 0, consistent with our findings in panel D) and d) of Fig.[9]for n > 2.5 galaxies. 
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With these results in mind, we now return to understand the changing location of 
^Mchar, the break in the SSFR— E M relation. Motivated by our SDSS analysis above, we 
show in Fig. [T2] the fraction of n < 1.5 and n > 2.5 Sersic index galaxies compared with 
the total number of galaxies (always with b/a > 0.55) in different Em bins for zCOSMOS 
and SDSS. The cross-over point Y*Mtrans shifts to higher surface mass densities, presumably 
because, while the average Em of low Sersic index galaxies (n < 1.5) remains unaltered, 
that of the higher Sersic index galaxies (n > 2.5) shifts at higher redshifts to higher values, 
because of the apparent size evolution noted above. This shift in T, Mtrans then produces the 
shift in YiMchar in the SSFR— Em relation. 

The behaviour of the SSFR - Em and stellar mass - size relations for 0.7 < z < 0.9 
zCOSMOS galaxies (Figs.E] and [TD]) in comparison with the equivalent SDSS samples (se- 
lected at higher masses because of the higher mass completeness limit in zCOSMOS, see 
Fig.H]) is similar to that of galaxies at the slightly lower redshift of 0.5 < z < 0.7, on which 
we have primarily focussed in the foregoing. However, the change in sizes of the early-type 
galaxies becomes even more pronounced at 0.7 < z < 0.9 and this leads to a further increase 
in Y^Mtrans arid E 

Mchar- 



4. Summary and conclusions 

We have studied the relationships between SSFR and stellar mass density in two mass- 
complete samples of galaxies (logM* > 10.4 for 0.5 < z < 0.7 and logM* > 10.7 for 0.7 < 
z < 0.9) drawn from the zCOSMOS survey, and the equivalent mass-complete samples of 
SDSS objects to draw the following conclusions: 



The median SSFR of SDSS galaxies is almost independent of Em for low values of 
Em, but then abruptly changes at a surface mass density Em c /w ~ 8.5, in agreement 
with previous studies. This step-function is clearly due to the change-over of different 
structural types from disk-dominated low Sersic galaxies (n < 1.5) to bulge-dominated 
high Sersic galaxies (n > 2.5) as the Em increases. The population mix changes over 
at a transition surface mass density logT>Mtrans ~ 8.45, which is almost identical to the 
characteristic sur f ace m ass density at which the SSFR changes, T,Mchar, identified by 



Kauffmann et all (120061 ) . 



The shape of the SSFR— Em relation in zCOSMOS at 0.5 < z < 0.9 is very similar to 
that of the SDSS, with a roughly uniform increase in the average SSFR by a factor of 
5-6 that is broadly independent of E M and which occurs in both early and late type 
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galaxies with n > 2.5 and n < 1.5 respectively. There is also a small increase of 0.1-0.2 
dex in T, Mchar at the higher redshifts. 

The rise in SSFR is almost exactly the same as that seen in the overall star-formation 
rate density (SFRD) of the Universe to this redshift, implying that galaxies across 
the full range of Em, arid with a wide range of Sersic index, are all contributing, 
proportionally, to the increase in the SFRD. 

The modest increase in TiMchar by ~ 0.1 — 0.2 dex in zCOSMOS relative to the SDSS 
is naturally explained by differences in the size-mass relations for disk-dominated (n < 
1.5) and bulge-dominated (n > 2.5) galaxies. Whereas the for mer have a si z e-mas s 



relation that does not change with redshift, as also found by iBarden et al.l (12 005). 



the la tter are smaller at higher redshifts, as also seen in the samples of iTrujillo et al. 



( 120071 ) . pushing individual galaxies to high E M values. This increases S Mtrans , and 



thus T.Mchar, at the higher redshifts. 

The median SSFR of disk-dominated galaxies (n < 1.5) is almost independent of 
surface mass density at both redshifts, but is about six times higher in zCOSMOS 
galaxies at z ~ 0.7 than in SDSS. With this strong decline in the average SSFR, the 
masses of the disk-dominated galaxies are unlikely to grow by more than ~ 30% from 
2; ~ 0.7 to 2 ~ 0. Coupled with the observed flat relation between size and stellar 
mass, this makes t he evidence for "ins ide-out" growth of disks weak, at least over this 



redshift range (c.f. IBarden et al.l 120051 ). 



The upturn in the median SSFR for n > 2.5 galaxies at lower Em is evidently due 
to galaxies with a significant disk component that have Sersic index n > 2.5 due to a 
large dominant bulge. Their exclusion results in a median SSFR that is independent 
of E^- also for the remaining n > 2.5 zCOSMOS galaxies, but is about a factor of five 
to six lower than for the n < 1.5 zCOSMOS galaxies. 



Putting these together we arrive at the following main observation concerning the evo- 
lution of relatively massive galaxies since z ~ 1. In parallel with the evidence that, at all 
redshifts, the mean SSFR within a given population (either disk-dominated with n < 1.5 
or bulge-dominated with n > 2.5) is independent of Em, and that the observed SSFR— Em 
step-function relation is due, at all redshifts, to the changing mix of disk-dominated and 
bulge-dominated galaxies as T> M increases and the strong difference in the average SSFR be- 
tween disks and bulges, we find that the increase in SSFR with redshift is also independent 
of Em and also of Sersic index n. The increase matches that of the global star-formation rate 
density of the Universe as a whole, suggesting that all types of galaxies are participating in 
the increase in SFR. 
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The conclusion is that the internal build-up of stellar mass in a galaxy is not strongly 
affected by E M , beyond the overall structural role of n (which clearly has a very large effect 
on the average SSFR, no doubt through the different structural components of disk and 
spheroid). The (relative) mean increase in SSFR back to z ~ 1 not only does not depend on 
Em, it does not even seem to depend on n. 
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Fig. 1.— Comparison of the stellar masses compute d using EqJT] with masses derived 
usi ng the entire COSMOS optical to infrared SED and iBruzual fc Charlotl (120031 ) models 
by iBolzonella et al. (boost ). A good agreement with a statistical rms of around 0.13 dex 
per galaxy and an offset (to higher SED masses) of 0.10 dex in the mean i s seen, reflecting 



typica l uncertainties in derived stellar masses. Interestingly, as described by Bolzonella et al. 



(120081 ). population synthesis models with TP-AGB phase (jMarastonl 120051 ) would produce 
a systmatic shift of ~ 0.1 dex towards lower SED masses than those shown in the figure, 
thus eliminating the offset between the masses computed using Eq.[T] and the SED masses. 
In order to achieve the highest possible internal consistency for our study, we derive mass 
estimates from purely rest-frame optical colors, for both SDSS and zCOSMOS, using Eq.[TJ 
(c.f. discussion in Sect. 12. 1.51) . 
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Fig. 2. — Comparison between stellar mass surface density and stellar mass for zCOS- 
MOS and SDSS mass-complete samples (c.f. FigHJ. The median E M values in different 
stellar mass bins are shown as filled circles, and 16th and 84th percentiles as solid lines. 
A slight increase of £ju with stellar mass, but with a large scatter, can be seen for both 
zCOSMOS and SDSS samples of relatively massive galaxies (c.f. Sect J2.1.5I) . 
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Fig. 3. — The stellar mass as a function of redshift for the SDSS sample (panels a and b), and for zCOSMOS 
galaxies at 0.5 < z < 0.9 (panels c-f) for different Sersic indices, n. The analysis in this paper is based on 
the mass-complete samples with logM* > 10.4 and logM* > 10.7 for 0.5 < z < 0.7 and 0.7 < z < 0.9. In 
order to be complete for SDSS galaxies (panel b) for the same mass range (logM* > 10.4) as for zCOSMOS 
galaxies, a redshift range of 0.04 < z < 0.08 for the comparison SDSS sample was chosen. 
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Fig. 4. — SSFR versus stellar mass as a function of Sersic indices for the 3048 zCOSMOS galaxies at 
0.5 < z < 0.9 selected as described in Sects. 12.1.11 and 12.1.21 The solid and dashed diagonal lines in each 
panel show star formation rates of 1 and 5 M Q /yr, respectively. The solid green lines show the paths of 
a model galaxy to redshift z ~ 0.6 (left panel) and z ~ 0.8 (right panel) with different redshifts Zf of the 
assumed onset of the star formation, and different e- folding timescale of the star formation t: Zf = 1 and 
r = 5Gyrs, Zf = 2 and r = 3Gyrs, Zf = 3 and r = 1 Gyr, and Zf = 5 and r = 1 Gyr, respectively 
(from left to right). The large magenta filled circles show the location of a set of models normalised to have 
Iab = 22.5 (the zCOSMOS se lection limit) at z = 0.7 (left panel), and z = 0.9 (right panel), produced 



Bruzual fc Charlotl \200j) models with a wide range of different star formation histories. These differ 



using 

in having different SFR c- folding timescales (r from 5 Gyrs to 0.7Gyrs), and different starting redshifts (zf 
from 1.3 to 5). The magenta filled circles effectively indicate the completeness limit of the sample, in the sense 
that the sample misses galaxies that lie to the left of these. The adopted mass completeness limits are shown 
by the green vertical thick lines at logM* > 10.4 for 0.5 < z < 0.7, and at logM„ > 10.7 for 0.7 < z < 0.9. 
This figure shows two well-known phenomena: first, the well-known downsizing effect (illustrated by the 
much higher number of galaxies in the yellow hatched region at 0.7 < z < 0.9 compared to 0.5 < z < 0.7), 
and secondly, the fact that galaxies with different Sersic indices have different star formation histories, even 
at a fixed mass (c.f. discussion in Sect l3.1j) . 



-30- 



t — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 

M ' 

x x 0.04<z<0.08 n<1.5 b/a>0.55 +>, 



logM>10.4 



zr o 



-l - 



-2 - 



-3 - 




zr o 



-l 



-2 - 



-3 L 




s 



N 



-H- 



+++ 



+4 



4- 



\ 

: . \ 



n — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 

N d) : 

x . 0.04<z<0.08 n>2.5 b/a>0.55 
N x logM>10.4 



\ 

\ 




- N b) 

N SDSS 0.04<z<0.08 logM>10.4 a x Q ; 04<z<0.08 b/a>0.55 logM>10.4 



X 



N 
X 

\ 



X 



X 




7 8 9 10 7 

log[£ M /[M kpc-2]] 



8 9 10 

log[S M /[M kpc-2]] 



Fig. 5. — Specific star formation rates for SDSS logM^ > 10.4 galaxies at 0.04 < z < 0.08 - the sample 
that should be directly comparable to the 0.5 < z < 0.7 zCOSMOS sample - are shown as green dots, the 
m edian SSFR values in di fferent Em bins as magenta triangles, and 25th and 95th percentiles (as in Fig. 9 
of Kauffmann et al. 200(j| ) as solid magenta lines. The diagonal dashed lines correspond to SFR surface 
densities ^sfr = 0.01Mo/yr/kpc 2 , and 0.1Mo/yr/kpc 2 , respectively. Panel a) shows all SDSS galaxies 
at 0.04 < z < 0.08 with logM„ > 10.4, while panel b) shows only the SDSS galaxies with an axis ratio 
b/a > 0.55 (see 12.1. 4]) . Panel c) shows the low Sersic index (n < 1.5) objects with b/a > 0.55, and panel d) 
the n > 2.5 SDSS galaxies, again with b/a > 0.55. The trend seen in panel a) and b) is clearly due to the 
change-over of different structural types from disk-dominated low Sersic index (n < 1.5) galaxies (in panel 
c) to bulge-dominated high Sersic (n > 2.5) index galaxies (panel d) as the Em increases. This change-over 
is also clearly seen in Fig. [51 
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Fig. 6. — The fraction of SDSS objects (with b/a > 0.55) with n < 1.5 (open magenta triangles) and with 
n > 2.5 (filled magenta triangles) in the logM + > 10.4 mass-complete sample. The fraction of n > 2.5 SDSS 
objects shows a sharp inc rease at logYiMtm n x ~ 8 .45, very close to the characteristic surface mass density 
logT, Mchar ~ 8.5 found bv lKauffmann et al.1 (|200q ). 
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Fig. 7. — The specific star formation rate SSFR versus stellar surface mass density Em of the zCOSMOS 
mass-complete sample at 0.5 < z < 0.7 (black dots), compared to the median and percentile curves sum- 
marizing the SDSS relations from Fig. [5] (magenta points and lines). The diagonal dashed lines correspond 
to SFR surface densities Ssfj?, = 0.01MQ/yr/kpc 2 , and 0.1MQ/yr/kpc 2 , respectively. The SSFR median 
values of zCOSMOS galaxies for different Em bins are shown as cyan filled squares, while for the highest Em 
we also show the mean SSFR values derived from co-added spectra as green open squares. Panel A) shows 
all galaxies in the zCOSMOS mass-complete sample, panel B) only galaxies with an axis ratio b/a > 0.55, 
panel C) galaxies with b/a > 0.55 and n < 1.5, and panel D) galaxies with b/a > 0.55 and n > 2.5. The 
respective small panels above the main ones (and the insert in panel A with finer grid in Em) show as filled 
squares (circles) the difference between the median zCOSMOS and SDSS SSFR values in a given Em bin. 
For n < 1.5 (disk) galaxies the SSFR stays roughly constant with Em, with the median SSFR at a given 
Em being about six times higher at z ~ 0.6 than in SDSS. Red arrows in panel c are model galaxy tracks 
discussed in Sect. 13.41 For n > 2.5 (early-type) galaxies the SSFR declines with Em for both zCOSMOS and 
SDSS galaxies, but there is again a shift to higher SSFR with about the same factor of about five to six, and 
also a modestly higher Em at the higher redshift. Almost all the galaxies causing the upturn at low Em 
are clearly galaxies with disks, and the median SSFR of n > 2.5 genuine early-type galaxies (see Sect. l3.3[) 
is also almost independent of Em (red squares in panel d) . 
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Fig. 8. — The specific star formation rate SSFR versus stellar surface mass density £m of the zCOSMOS 
mass-complete sample at 0.7 < z < 0.9 (black dots), compared to SDSS (magenta symbols and lines). 
Symbols are as in Fig. [71 and similar trends of the SSFR as a function of Em and Sersic index for 0.7 < z < 0.9 
zCOSMOS galaxies are seen as for the slightly lower redshift 0.5 < z < 0.7 objects (Fig. [7]). 
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Fig. 9. — The semi-major half-light radius r 1 / 2 in the observed ACS I-band (corresponding roughly to 
rest-frame B-band for the redshift range 0.5 < z < 0.9) versus stellar mass for the zCOSMOS mass-complete 
sample at 0.5 < z < 0.7. Individual zCOSMOS measurements are shown as black dots, and median values as 
filled cyan squares. The median SDSS half-light radii in different mass bins derived from g-band images are 
shown as filled magenta triangles, while the 16th and 84th percentiles of r 1 / 2 of the SDSS and zCOSMOS 
galaxies in each mass bin are shown as solid magenta and cyan lines respectively. The two dashed diagonal 
black lines show mass surface densities of log^M = 9 (upper line), and log^M = 9.5, respectively. Panel a) 
shows all galaxies in the mass-complete sample, panel b) galaxies with an axis ratio b/a > 0.55, panel c) 
galaxies with b/a > 0.55 and n < 1.5, and panel d) galaxies with b/a > 0.55 and n > 2.5. The respective 
small panels above the main ones show as filled squares the difference between the median zCOSMOS and 
SDSS Ti/2 values in a given mass bin. For n < 1.5 galaxies there is almost no evolution in the stellar mass - 
size relation between z ~ 0.7 to z ~ 0, while for n > 2.5 objects the average half-light radius of galaxies at 
a given mass is smaller by ~ 25% at z ~ 0.7. 
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Fig. 10. — The semi-major half-light radius rw2 versus stellar mass for the zCOSMOS mass-complete 
sample at 0.7 < z < 0.9 (black dots and cyan symbols and lines) compared to SDSS (magenta symbols 
and lines). Symbols are as in Fig. [5J and similar trends of the stellar mass - size relation are seen here for 
0.7 < z < 0.9 zCOSMOS galaxies as seen in Fig.[9] for the slightly lower redshift 0.5 < z < 0.7 objects. 
However, the change in sizes of the early-type galaxies (panel d) becomes even more pronounced and this 
leads to a further increase in trans and Tim char- 
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Fig. 11. — Comparison of the B-band M/L ratios for zCOSMOS galaxies with b/a > 0.55 and n > 2.5 
(black filled circles) w ith the M/L ratio s obtained studying the fundamental plane for early-type (E+SO) 
galaxies from Fig. 15 of Treu et al. ( 2005) , and the local relation for early-type (spheroid) galaxies taken from 
the same figure of iTreu et al.l (|2005l dashed line). The slope of the observed mean zCOSMOS M/L-mass 
relation agrees with the slope of the local relation, as indicated by the solid diago nal line, offset to l ower M/L, 
as expected for a passively evolving population. Moreover, the spheroids from I Treu et alJ ( 2005h occupy a 
similar region of the M/L-mass diagram as the zCOSMOS galaxies at 0.5 < z < 0.7. This reassures that the 
derived stellar masses for n > 2.5 zCOSMOS galaxies are reasonable, and further suggests that the changes 
in the stellar mass-size relation at 0.5 < z < 0.7 in Fig. [9] are due mainly to smaller sizes (or at least smaller 
measured rt/2), at a given mass, in zCOSMOS compared with SDSS. 
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Fig. 12. — The fraction of b/a > 0.55 objects with n < 1.5 (blue filled squares) and n > 2.5 (red filled 
squares) for zCOSMOS galaxies at 0.5 < z < 0.7 from the mass-complete sample with logM„ > 10.4. The 
respective fractions for SDSS galaxies are shown as filled and open magenta triangles, for n > 2.5 and n < 1.5 
galaxies, respectively. The transition surface mass density TtMtrans (dot-dashed cyan vertical line) at which 
these curves cross is larger for zCOSMOS than for SDSS galaxies (dashed magenta vertical line), and gives 
an explanation for the shift seen in the location of the step "Sivichar in the median SSFR versus Em relation 
for zCOSMOS galaxies compared to that of the SDSS (Figs. [7] and [8]). 
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A. zCOSMOS vs. SDSS size measurements 



It might be a concern that there are significant systematic differences between the size 
measurement methods for zCOSMOS (as described in Sect. 12. 1.31) and SDSS (as described in 
Sect. 12.2b ). To check the consistency of our size measurements for the two samples we used 
the ACS images of the mass-complete sample of ~ 1200 zCOSMOS galaxies at 0.5 < z < 0.9 
to test two issues: i) are the sizes determined using a circular aperture (as in the case of SDSS) 
with a correction of yajb consistent with the sizes determined using elliptical apertures (as 
in GIM2D applied to zCOSMOS); and ii) are the GIM2D half-light radii consistent with 
a derivation of half-light radii starting from the Petrosian radius containing 50% of the 
Petrosian flux and applying the corrections done for the SDSS sample described in Sect J2. 2a 
? The results are shown in the two panels of Fig. [131 an d show that there are no significant 
systematic differences between the two size measurement methods in zCOSMOS and SDSS. 
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Fig. 13. — The left panel shows a good agreement (offset 0.009 dex, RMS 0.03 dex) between 
the COSMOS sizes (semi-major axis) computed directly from the ACS images using an 
elliptical aperture, and the sizes computed using circular apertures and multiplying them by 
a/ a/b. The right panel shows also a quite good agreement (offset 0.008 dex, RMS 0.07 dex) 
between the sizes (semi-major axis) computed with GIM2D, and the sizes computed starting 
fro m the circular Petrosia n radius containing 50% of the Petrosian flux and applying equation 
6 of lGraham et al.l ( 20051 ) and the y/ a/b correction to compute the semi-major half-light radii 
(the same method as applied to SDSS and described in Sect. 12. 21 a). 



